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Abstract X-ray binaries are excellent laboratories to study collapsed objects. On 
the one hand, transient X-ray binaries contain the best examples of stellar-mass 
black holes while persistent X-ray binaries mostly harbour accreting neutron stars. 
The determination of stellar masses in persistent X-ray binaries is usually hampered 
by the overwhelming luminosity of the X-ray heated accretion disc. However, the 
discovery of high-excitation emission lines from the irradiated companion star has 
opened new routes in the study of compact objects. This paper presents novel tech- 
niques which exploits these irradiated lines and summarises the dynamical masses 
obtained for the two populations of collapsed stars: neutron stars and black holes. 



1 Introduction 

X-ray binaries (XRBs hereafter) are interacting binaries where X-rays arise from the 
accretion of matter onto a neutron star (NS) or a black hole (BH). Accretion pro- 
cesses are found in other astrophysical environments such as cataclysmic variables 
(i.e. interacting binaries with accreting white dwarfs), T Tauri stars, protoplanetary 
discs, etc. but the unique property of XRBs is the presence of central compact ob- 
jects that are the remnants of collapsed, massive stars. Therefore, they provide the 
best laboratories to study their properties in detail, such as masses, spin or NS equa- 
tion of state. This paper is not meant to give a thorough review of XRBs but instead 
it will focus on three selected topics with implications for our knowledge of the 
mass spectrum of collapsed stars. These are: 

1 . Evidence for BHs in XRBs: a summary of dynamical masses is presented. 

2. The Bowen Project: a new technique to trace the orbit of companion (donor) 
stars in persistent XRBs. 
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3. Echo Tomography: reprocessed light from the donor star can constrain the bi- 
nary inclination, the parameter with largest impact in the mass error. 

Excellent reviews on other aspects of XRBs can be found in 191 and ll24l . 



2 Black Holes in X-ray Binaries 

The mass distribution of Black Holes (BHs hereafter) has strong impact in sev- 
eral areas of Astrophysics, in particular SNe models, the evolution of massive stars, 
chemical enrichment of the galaxy, jet formation etc.. Stellar evolution theories pre- 
dict ~ 10 9 BHs remnants in the Galaxy (2) but only BHs in interacting binaries can 
be easily detected through X-ray radiation triggered by accretion. This is the reason 
why the history of BH discoveries has run in parallel with the development of X-ray 
astronomy. 
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Fig. 1 Cumulative distribution of 
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In particular Soft X-ray transients (SXTs hereafter) provide the best systems to 
find stellar-mass BHs, since ~75% of these transients likely harbour a BH. SXTs 
are a subclass of X-ray binaries with low-mass donor stars (typically of K-M spec- 
tral types) which exhibit episodic outbursts due to mass transfer instabilities in the 
accretion disc |23l . During outburst, the X-ray luminosity increases by factors up 
to 10 6 -10 7 and, therefore, they are easily spotted by X-ray satellites. Unfortunately, 
the companion star is overwhelmed by the X-ray heated disc at all wavelengths, 
precluding its detection. However, after a few months of activity, X-rays switch off, 
the reprocessed light drops off several magnitudes into quiescence and the com- 
panion starts to dominate the optical flux. This provides a special opportunity for 
spectroscopic detection, perform dynamical studies and derive stellar masses. Fig- 
ure 1 shows the cumulative histogram of BH SXTs discovered since 1966, when 
Cen X-2 was first detected during a rocket flight. A linear increase at a rate of ~ 1.7 
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yr is apparent since the late 80's, when a new fleet of X-ray satellites with higher 
sensitivity and All-Sky-Monitor capabilities became operative. 

The best evidence for the presence of BHs is dynamical, i.e. a compact object 
whose mass exceeds the maximum allowed for stable neutron stars or ~3 M l36l . 
And this is relatively easy to prove through the observed radial velocity curve of 
the companion star during quiescence. The orbital period P or b and velocity semi- 
amplitude Kq combine in the standard mass function equation f(M) 

f(M\ - Po * K c - M *sin 3 / 

which relates the mass of the compact object M\ and the companion starMc through 
the orbital inclination i and mass ratio q — ^ . It is easy to show that f(M) yields a 
secure lower limit to the BH mass Mx- This experiment typically requires resolving 
powers > 1500 and can be performed with targets brighter than R ~23 using current 
instruments on lOm-class telescopes. 

But getting actual BH masses rather than lower limits is not so straightforward. 
The reason being that, by their very nature, BHs do not burst nor pulse and hence one 
cannot trace their orbital motion. We are facing a single-line spectroscopic binary 
where the extra observables (q and i) must be extracted from the optical star. And 
this can be accomplished with two further experiments: 

• Resolving the rotational broadening of the donor's photospheric lines. Because 
the donor star is filling the Roche lobe and synchonized, the rotational broadening 
V sin i scales with Kq as a function of q ll42l . Therefore, by measuring V sin i one 
can determine q directly. 

• Fitting synthetic models to the ellipsoidal modulation. The changing visibility of 
the tidally distorted companion star generates a double-humped light curve (the 
so-called ellipsoidal modulation) whose amplitude is a strong function of i and 
q. For extreme mass ratios q < 0.2, typical of BH binaries, the shape of the light 
curve is weakly sensitive to q and hence ; can be easily determined [38]. 

By combining the mass function with constraints on q and i one gets a full dy- 
namical solution and hence the BH mass with minimum assumptions. Further de- 
tails on this prescription and possible systematics involved can be found in several 
reviews such as 0. 

BHs have also been found in a handful of High Mass X-ray Binaries (HMXBs 
hereafter), i.e. XRBs with early-type massive donor stars. However, here we find 
several limitations which complicate the analysis. A key factor is Mq, which for a 
HMXB is a large number and has a wide range of uncertainty. The optical star is 
likely to be undermassive for its spectral type as a result of mass transfer and binary 
evolution 1351 . Furthermore, mass transfer is usually produced through winds rather 
than Roche lobe overflow and this has a two-fold effect. On one side, wind emission 
can contaminate the radial velocities of the donor star. On the other, since the optical 
star does not fill its Roche lobe, q and i values derived through V sin / and ellipsoidal 
fits may be overestimated. These caveats can only be side-lined in eclipsing binaries 
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such as X-7 in M33, which so far is the only case where this has been possible. In 
addition to the eclipse duration, the distance provides an extra restriction which lead 
to tight constraints in the space parameter. In particular, the radius of the donor, the 
Roche lobe filling factor and the inclination are accurately determined and yield a 
very precisse BH mass [33). 



Table 1 Dynamical BHs 



System 


forb 


f( M ) 


Donor 


Classification 


Af x + 




[days] 


[M Q ] 


Spect. Type 




[M s ] 


GRS 1915+105° 


33.5 


9.5 ± 3.0 


K/M III 


LMXB /Transient 


14 ± 4 


V404 Cyg 


6.471 


6.09 ± 0.04 


K0IV 




12 ±2 


Cyg X-l 


5.600 


0.244 ± 0.005 


09.7 lab 


HMXB/Persistent 


10 ±3 


LMC X-l 6 


3.909 


0.143 ±0.007 


07 III 




10.3 ± 1.3 


M33 X-7 


3.453 


0.46 ± 0.08 


07-8 III 




15.7 ± 1.5 


XTEJ1819-254 


2.816 


3.13 ±0.13 


B9III 


IMXB/Transient 


7.1 ± 0.3 


GRO J1655-40 


2.620 


2.73 ± 0.09 


F3/5 IV 




6.3 ± 0.3 


BW Cii x 


2.545 


5.73 ± 0.29 


G5IV 


LMXB/Transient 


> 7.0 


GX 339-4 rf 


1.754 


5.8 ±0.5 






> 6.0 


LMC X-3 


1.704 


2.3 ±0.3 


B3 V 


HMXB/Persistent 


7.6 ± 1.3 


XTEJ1550-564 


1.542 


6.86 ± 0.71 


G8/K8 IV 


LMXB/Transient 


9.6 ± 1.2 


4U 1543-475 


1.125 


0.25 ± 0.01 


A2 V 


IMXB/Transient 


9.4 ± 1.0 


H1705-250 


0.520 


4.86 ±0.13 


K3/7 V 


LMXB/Transient 


6±2 


GS 1124-684 


0.433 


3.01 ±0.15 


K3/5 V 




7.0 ± 0.6 


XTE J 1859+226'" 


0.382 


7.4 ± 1.1 








GS2000+250 


0.345 


5.01 ±0.12 


K3/7 V 




7.5 ± 0.3 


A0620-003 


0.325 


2.72 ± 0.06 


K4V 




11 ±2 


XTE J 1650-500 


0.321 


2.73 ± 0.56 


K4V 






GRS 1009-45 


0.283 


3.17 ± 0.12 


K7/M0 V 




5.2 ± 0.6 


GRO J0422+32 


0.212 


1.19 ±0.02 


M2 V 




4± 1 


XTE Jill 8+480 


0.171 


6.3 ± 0.2 


K5/M0 V 




6.8 ± 0.4 



' Masses compiled by (9) and |32| . 

" New photometric period of 30.8±0.2 days reported by 1301 . 

* Updated after (34). 

Updated after (7). 

d Updated after (27). 

" Period is uncertain. See 1431 . 

Table Q] presents an updated list of confirmed BHs based on dynamical argu- 
ments, with their best mass estimates. We currently have 21 BHs, with orbital pe- 
riods between 33.5 days and 4.1 hours. The great majority are SXTs (17) while 4 
are persistent HMXBs: Cyg X-l plus the 3 extragalactic binaries LMC X-l, LMC 
X-3 and M33 X-7. The case of GX 339-4 deserves special mention because it is the 
only SXT where the presence of a BH was proven during the outburst phase. This 
was possible thanks to the detection of fluorescent lines arising from the irradiated 
companion (see Sect. 01. GRS 1915+105 is also noteworthy, not only because of its 
long orbital period and large mass function but also because IR spectroscopy was 
essential to overcome the > 25 magnitudes of optical extinction and reveal the radial 
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velocity curve of the companion star |20l . However, it should be noted that recent 
photometry reports a slightly shorter orbital period and evidence for irradiated light 
curves 11301 . The combination of these two effects will likely decrease the mass 
function and BH mass. XTE J 1859+226 also needs revisiting because its orbital pe- 
riod is uncertain J43). In summary we have 16 BH masses ranging between 4 and 
16 M© with ~5-30 % errors. These can be compared with theoretical distributions 
of stellar remnants such as 1191 . The model includes binary interaction under Case 
C mass transfer (i.e. Common Envelope evolution after core helium ignition), wind 
mass-loss in the Wolf-Rayet phase and SN lb explosion. The computation predicts 
a continuum distribution of remnants with a mass cut at 12 M© which is difficult to 
reconcile with some of the observed masses. However, the model entails many the- 
oretical uncertainties which dominate the final mass spectrum such as the Common 
Envelope efficiency, the wind mass-loss rate or the progenitor's mass cut. Clearly 
more SXT discoveries and lower uncertainties in BH masses are required before 
these issues can be addressed and the form of the distribution is used to constrain 
BH formation models and XRB evolution. 
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Fig. 2 Magnitude distribution of 
BH SXTs in quiescence. The black 
histogram indicates dynamical BHs 
while the rest are BH candidates. 
Targets without optical counterpart 
are likely fainter than R > 23. 



In addition to dynamical BHs, there are 27 other SXTs with similar X-ray spec- 
tral and timing properties during outburst Unfortunately, these BH candidates be- 
come too faint in quiescence for dynamical studies or even lack accurate astrometry. 
This is illustrated in Fig. 2 which shows the magnitude distribution of the 44 cur- 
rently known BH transients. Dynamical studies are only possible with the current 
largest telescopes for sources brighter than R < 23. Not shown in the figure is the 
heavily reddened GRO 1915+105 which was studied in the NIR. The figure de- 
picts the bright tail of a dormant population of galactic BH SXTs which several 
works have estimated in a few thousand systems ( |3?! 1 and included references). Im- 

1 This number has been updated after 1241 with new detections reported in several Astronomical 
Telegrams. 
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proving the statistics of dynamical BHs requires not only a new generation of ELT 
telescopes to tackle fainter targets but also new strategies aimed at unveiling new 
"hibernant" SXTs before they go into outburst. Quiescent BH SXts typically have 
EW (H a )~ 20 — 50 A and hence they should show up in deep H a surveys such as 
IPHAS ifTTl . However, clever diagnostics need to be defined to clear out other popu- 
lations of H a emitters such as cataclysmic variables or T Tauris (see Corral-Santana 
et al., these proceedings). 



3 The Bowen Project 

Aside from transient XRBs, there are ~150 persistent XRBs in the Galaxy, the 
great majority hosting neutron stars (NS hereafter) accreting at the Eddington limit. 
They are considered the progenitors of Binary Millisecond Pulsars (BMPs hereafter) 
because is the sustained accretion during their long active lives that spins the NS up 
to millisecond periods. The discovery of millisecond pulses in 8 transient XRBs and 
coherent oscillations during X-ray bursts in 13 persistent XRBs gave strong support 
to this "recycle" pulsar scenario. And burst oscillations were detected in addition 
to persistent pulses in the transient XRBs SAX J1808-3658 (TO] and XTE J1814- 
338 l40l with identical frequencies. This confirmed that burst oscillations are indeed 
modulated with the spin of the NS. The interest of these discoveries stands in the 
fact that one can use the orbital Doppler shift of pulses/oscillations to trace the NS 
orbit and obtain the X-ray mass function. 

Optical emission in persistent XRBs is triggered by reprocessing of the intense 
X-ray radiation in different binary sites, mainly the accretion disc. The compan- 
ion star is ^1000 times fainter than the irradiated disc at optical-IR wavelengths 
and hence completely undetectable. This has systematically plagued attempts to de- 
termine system parameters and, in most cases, only the orbital period is known. 
Fortunately, there are methods which can exploit the effects of irradiation and X-ray 
variability. New prospects were opened by the discovery of sharp high excitation 
emission lines arising from the irradiated face of the companion star in Sco X- 1 ||39ll . 
The most prominent are found in the core of the Bowen feature, a blend of CIII/NIII 
lines which are mainly powered by fluorescence. These lines trace the motion of 
the companion star and provided the first dynamical information on this protypical 
LMXB (see Fig. 3). Since then, sharp Bowen lines from companion stars have been 
discovered in 7 other persistent LMXBs and 4 transients during outburst: Aql X-l, 
GX 339-4 and the BMPs XTE J1814-338 and SAX J1808.4-3658. These transient 
studies beautifully demonstrate the power of this technique in systems which oth- 
erwise cannot be studied in quiescence because either they are too faint (case of 
GX 339-4 and the BMPs) or are contaminated by a bright interloper (Aql X-l). In 
particular, the case of GX 339-4 is remarkable because the Bowen study provides 
the first solid evidence for the presence of a BH in this classic transient. 

The radial velocity curves of the Bowen lines are biased because they arise from 
the irradiated face of the star instead of its center of mass. Therefore, a K-correction 
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Fig. 3 Detecting companion stars in persistent XRBs. Left: the main high excitation emission lines 
due to irradiation of the donor star in Sco X-l. Adapted from 1391 . Right: Trail spectrum showing 
the radial velocity motion of the Bowen CIII/NIII lines as a function of time. After 1391 . 

needs to be applied in order to obtain the true velocity semi-amplitude Kq from 
the observed velocity K em . The K-correction parametrizes the displacement of the 
center of light with respect to the donor's center of mass through the mass ratio and 
disc flaring angle a. The latter dictates the size of the disc shadow projected over the 
irradiated donor [25). Extra information on q and a is thus required to get the real 
Kq ■ Furthermore, useful limits to the NS mass can be set if the binary inclination is 
well constrained through eclipses. 

Table [2] summarises the NS masses obtained through the Bowen technique dur- 
ing several campaigns at the WHT, AAT and VLT. The list of persistent systems is 
almost a complete sample of Galactic LMXBs brigther than B~19. In the cases of 
Aql X-l and XI 822-371 the evidence of NS more massive than canonical is very 
persuasive. The latter is a particularly favourable binary because it is eclipsing and 
the NS is a pulsar. Then its radial velocity curve is known through the study of or- 
bital pulse delays. Good constraints on the NS velocity are also available for V801 
Ara through the detection of pulse oscillations during a superburst J6). Tight limits 
to the inclination and mass ratio are also available for the eclipsing EXO 0748- 
676 l29l and the dipper GR Mus 0"). In the remaining cases the NS mass is not well 
constrained due to large uncertainties in the inclination and/or mass ratio. However, 
it is important to stress that these are the first dynamical constraints in persistent 
LMXBs since their discovery, 40 years ago. Other techniques (such as the Echo 
Tomography) need to be exploited to further refine these limits and derive more ac- 
curate NS masses. Previous reviews presenting results of the Bowen project can be 
found in 151 and fl4l. 
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Table 2 NS masses obtained using the Bowen Technique 



System 


Porb 


Mag 


Type 






Reference 




[hr] 






[km/s] 


[M Q ] 




Sco X-l 


18.9 


B=12.2 


Persistent 


87 ± 1 


>0.2 


(39) 


LMC X-2 


8.1 


B=18 




351 ±28 


> 1.2 


fT3l 


X1822-371 


5.6 


B=15.8 




300 ± 15 


1.6-2.3 


ID 


V926 Sco (X1735-444) 


4.7 


B=17.9 




226 ± 22 


> 0.5 


m 


GX9+9 (X1728-16) 


4.2 


B=16.8 




230 ± 35 


> 0.3 




GR Mus 


3.9 


B=19.1 




245 ± 30 


1.2-2.6 


m 


V801 Ara(X1636-536) 


3.8 


B=18.2 




277 ± 22 


> 0.8 





EXO 0748-676 


3.8 


B=16.9 




310 ± 10 


1.1-2.6 


(29) 


Aql X-l 


19 


V-22 


Transient 


247 ±8 


>1.6 


CD 


GX 339-4 


42.1 


V >21 




317 ± 10 


>6.0 


ED [27] 


XTE J1814-338" 


4.2 


V= 23.3 


Transient/BMP 


345 ± 19 


>1.0 


181 


SAX J 1808. 4-365 8 


2.0 


V=21.8 




248 ± 20 


>0.3 


(2) 



1 After K-correction and constraints to the inclination, mass ratio or NS velocity (when available). 
" Preliminary results. 



4 Echo Tomography 

Echo Tomography uses time delays between X-ray and UV/optical variability as a 
function of orbital phase to map the reprocessing sites in a binary RD . The optical 
variability can be modelled by the convolution of the X-ray light curve with a trans- 
fer function which depends on the binary geometry. The transfer function encodes 
information on the most fundamental parameters such as the binary inclination, star 
separation and mass ratio. And in particular, the component associated with the 
companion star is most sensitive to these parameters so detecting echoed emission 
from the donor offers the best opportunity to constrain them. There has been sev- 
eral attempts at detecting correlated optical and X-ray variability using white light 
or broad band filters (e.g. BTI . l22l0 . These works have detected delays which are 
mostly consistent with reprocessing in the outer disc implying that the disc is the 
dominant source of continuum reprocessed light. 

Exploiting emission-line reprocessing rather than broad-band photometry has 
two potential benefits: a) it amplifies the response of the donor's contribution by 
suppressing most of the background continuum light (dominated by the disc); b) 
since the emission line reprocessing time is instantaneous, the response is sharper 
(i.e. only smeared by geometry). Through the Bowen project we know that high en- 
ergy radiation is very efficiently reprocessed by the donor's atmospheres into Bowen 
fluorescence lines. Therefore, we decided to search for optical echoes of X-ray vari- 
ability using ULTRACAM lfl6ll equipped with a special set of narrow band filters, 
centered at the Bowen blend and a red continuum. The latter is essential to subtract 
the continuum light and hence amplify the reprocessed signal from the companion. 
During an RXTE/WHT campaign on Sco X-l correlated variability was detected at 
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Fig. 4 Echo Tomography experiment in Sco X-l. Left: large amplitud X-ray variability and cor- 
related optical light curve observed at orbital phase 0.5. Sco X-l happened to be in the flaring 
branch state. Right: Cross-correlation functions between X-rays and optical light curves observed 
in the continuum (top), Bowen+Hell window (middle) and Bowen+Hell after continnum subtrac- 
tion (bottom). After (26). 



phase ~ 0.5 i.e. superior conjunction of the companion star, when the heated face 
presents its maximum visibility J26). Time delays of 14-16s are measured after the 
continuum light is subtracted from the Bowen light curves (see Fig. 4). These de- 
lays are consistent with the light traveltime between the NS and the companion star 
and hence provide the first evidence of reprocessing in the companion of Sco X-l. 
However, one needs to detect several optical echoes as a function of orbital phase in 
order to constrain i and q and derive masses. 

In a second campaign we observed the burster X1636-536 simultaneously with 
RXTE and VLT+ULTRACAM. Three X-ray bursts and their corresponding opti- 
cal echoes were recorded at orbital phases 0.55, 0.20 and 0.83 and these are shown 
in the left panel of Fig. 5. The optical bursts clearly lag X-ray burst and are also 
smeared, indicating an extended reprocessing site. Delay times are in the range 2-3 
sees showing little evidence for orbital variability. However, these delays drift when 
several amounts of continuum light (parametrised by the factor cf) are subtracted 
from the Bowen+Hell light. And for cf ~ 0.8 — 0.95 the 3 delays become consis- 
tent with reprocessing in the companion for Mns = 1.4 M Q , q = 0.3, a = 12° and 
i = 36 — 60° , as derived through radial velocities of the Bowen lines |6|. This is 
illustrated in the right panel of Fig. 5. Note that, in particular, delays observed at 
phase ~ 0.5 are especially sensitive to the inclination angle. The main difficulty 
which hinders us from constraining the inclination is the unknown amount of con- 
tinuum substraction. In principle, there must be an optimum cf factor which results 
in a perfect subtraction. However, this is not easy to find because the continuum 
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Fig. 5 Echo Tomography of X1636-536. Left: the three X-ray bursts detected and their optical 
(Bowen+Hell) counterparts. Right: delay times between the X-ray and Bowen+Hell burst light 
curves as a function of continuum subtraction factor. Shaded regions correspond to delays expected 
for reprocessing in the companion at each orbital phase. They are computed for M^s = 1.4 M , 
q = 0.3, i = 36 - 60° and a = 12°. After |28). 



filter is placed ~1500 A away from the Bowen lines due to the optical layout of 
ULTRACAM. New high-speed spectrophotometry devices such as ULTRASPEC 
will provide pure emission line light curves for echo mapping experiments. These 
are likely to yield accurate inclinations and, when combined with dynamical infor- 
mation from the Bowen lines and X-ray mass functions, the first accurate NS masses 
in persistent XRBs. 



5 Conclusions 



In the past 20 years the field of X-ray binaries has experienced significant 
progress with the discovery of 17 new BHs and 8 transient BMPs in LMXBs. 
Dynamical masses are available for 16 BHs but better statistics and improved 
errors are required before using the observed distribution to constrain XRB 
evolution and supernova models. Exploiting deep Ha surveys of the Galactic 
plane, such as IPHAS, may unveil a significant fraction of a large expected 
population of quiescent XRBs. 
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The discovery of fluorescence emission from the companion star has 
opened the door to derive NS masses in persistent and new transient XRBs. 
This is possible thanks to: i) dynamical information from irradiated donors 
through high-resolution spectroscopy of the Bowen blend; ii) echo-mapping 
reprocessing sites through simultaneous Bowen-line/X-ray lightcurves. These 
techniques, together with results from burst oscillations and transient BMPs, 
will likely provide the first accurate NS masses in XRBs in the near future 
and perhaps confirm the existence of massive NS. Thanks to these new tech- 
niques, which have proven their worth, the future is bright as new instruments 
and telescopes will allow to push ahead our sample of BHs and NS masses. 
High-speed and high-resolution instruments, such as OSIRIS at GTC, RSS at 
SALT and ULTRASPEC, will play a crucial role in this goal. 
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